Fate of oscillating scalar fields in the thermal bath 
and their cosmological implications 

Jun'ichi Yokoyama 

Department of Earth and Space Science, Graduate School of Science, 

Osaka University, Toyonaka 560-0043, Japan 

(Dated: February 2, 2008) 

O ; Abstract 

^^ ' Relaxation process of a coherent scalar field oscillation in the thermal bath is investigated using 

^ ■ nonequilibrium quantum field theory. The Langevin-type equation of motion is obtained which 

C^^ . has a memory term and both additive and multiplicative noise terms. The dissipation rate of 

the oscillating scalar field is calculated for various interactions such as Yukawa coupling, three- 
body scalar interaction, and biquadratic interaction. When the background temperature is larger 
than the oscillation frequency, the dissipation rate arising from the interactions with fermions is 



> 

(N 

o 

T-|- ! suppressed due to the Pauli blocking, while it is enhanced for interactions with bosons due to the 

induced effect. In both cases, we find that the microphysical detailed balance relation drives the 
oscillating field to a thermal equilibrium state. That is, for low-momentum modes, the classical 
O , fluctuation-dissipation theorem holds and they relax to a state the equipartition law is satisfied, 

while higher-momentum modes reach the state the number density of each quanta consists of 



p. 



c^ 



/\^ • the thermal boson distribution function and zero-point vacuum contribution. The temperature- 

dependent dissipation rates obtained here are applied to the late reheating phase of inflationary 
universe. It is found that in some cases the reheat temperature may take somewhat different value 
from the conventional estimates, and in an extreme case the infiaton can dissipate its energy without 
linear interactions that leads to its decay. Furthermore the evaporation rate of the Affleck-Dine 
field at the onset of its oscillation is calculated. 

PACS numbers: 98.80.Cq,11.10.Wx,05.40.-a OU-TAP-232 



I. INTRODUCTION 

Cosmology of the early Universe is a useful probe of high energy phenomena beyond the 
reach of ground based accelerator experiments. The universe at its birth, however, is likely 
to suffer from huge relic quantum fluctuations and we cannot expect that it started classical 
evolution from a thermal equilibrium state with a well-defined temperature. Rapid cosmic 
expansion in the early universe further delays equilibration [jj and it is not likely that the 
phase transition of grand unified theories occurred thermally [2]. Once the energy scale 
has fallen well below typical grand unification scale, cosmic expansion rate gets smaller than 
interaction rates of ambient massless particles to establish thermal equilibrium. Phenomenon 
in such a regime may be studied in terms of quantum field theory at finite temperature 
neglecting cosmic expansion and using the cosmic temperature at each epoch. If some 
degrees of freedom are out of equilibrium, then we must of course use nonequilibrium field 
theories |3|. In modern cosmology, we often encounter a situation some scalar fields are in 
nonequilibrium configuration interacting with thermal background. 

Indeed scalar fields play central roles to explain virtually everything we observe — overall 
homogeneity and isotropy as well as the origin of small density perturbation are attributed 
to inflation driven by an inflaton scalar field, huge entropy carried by the cosmic microwave 
and neutrino) background radiation to the reheating process by the decay of the inflaton 
m Furthermore the observed baryon asymmetry a,rd dark matter may also or,gmate m 
scalar fields such as squarks and/or sleptons through the Affleck-Dine mechanism jy] and 
formation of Q-balls [7 1 . 

Thus it is of utmost importance to clarify the evolution of scalar fields in cosmic medium. 
In the present paper we study the fate of a coherent scalar field oscillation interacting with 
fermions or bosons, which are thermally populated, using the nonequilibrium quantum field 
theory. Such a situation is realized in the late stage of reheating after inflation as well as 
in the evolution of flat directions in supersymmetric theories which may be associated with 
Affleck-Dine baryogenesis. 

We start with a brief review of a fleld theoretic method appropriate to analyze time 
evolution of the expectation value of a scalar fleld. The standard quantum fleld theory, 
which is appropriate for evaluating the transition amplitude from an 'in' state to an 'out' 
state for some fleld operator O, (out| O |in), is not suitable to trace time evolution of an 



expectation value in a non-equilibrium system. In order to follow the time development of 
the expectation value of some fields, it is necessary to establish an appropriate extension 
of the quantum field theory, which is often called the in-in formalism. This was first done 
by Schwinger p and developed in J9l [lO, lul- This method has been applied to various 
cosmological problems by a number of authors jl2, ll3|, [ij, ll5|, lla llTl llSl Il9l |20|, |2l| . To 



name a few, a Langevin equation has been obtained by Morikawa [1^ and Gleiser and Ramos 
[15[ in the slow- roll limit, which was applied to the electroweak phase transition in |l7| and 
to warm inflation 22 1 in |20|. On the other hand, the case of oscillating scalar field was 

n 

studied by Greiner and Miiller who took only the self interaction into account [1SI|. Our 



work is partially related to it but we consider more general interactions with other fermions 
and bosons, whose effects are strikingly different from each other as shown in [l7[. 

We calculate an effective action for a real scalar field (p perturbatively in the in-in formal- 
ism by integrating out fields interacting with assuming that they are in thermal equilibrium 
distributions at a fixed temperature in a fixed flat spacetime. The resultant effective action 
is complex-valued as a result of coarse graining of these interacting fields, and it describes 
dissipation of the system field 0. This complex-valuedness is cured by the introduction of 
auxiliary fields which act as noise terms, both additive and multiplicative, in the equation 
of motion. Its derivation from the effective action is reported in the next section. 

In §111 the equation of motion is explicitly solved in the case only linear terms in </> 
are important. We show that each spatial Fourier mode of the scalar field will relax to 
a value determined by the ratio of the Fourier transform of the noise correlation function 
and that of the memory kernel in the equation of motion, and it takes the same thermal 
equilibrium value for all the three interactions discussed there, namely Yukawa coupling, 
three-body scalar interaction, and biquadratic interaction. This is achieved by the detailed 
balance relation which also leads to the classical fluctuation-dissipation theorem for low 
momentum modes. The time scale for the relaxation, which is essentially important for 
cosmological applications, is also evaluated for respective interactions. The result is quite 
different depending on the statistical property of the interacting particles. 

In §IV the analysis is extended to the multiplicative noises and dissipation. Although 
we cannot find a solution to the equation of motion in this case, we can still confirm the 
generalized fluctuation-dissipation relation and obtain the dissipation rate as well. 

These formulae are applied to two cosmological situations, namely the late reheating 



phase after inflation in §V and oscillating flat direction in §VI. Finally §VII is devoted to 
summary and discussion. 



II. EFFECTIVE ACTION IN NONEQUILIBRIUM QUANTUM FIELD THEORY 

A. Nonequilibrium quantum field theory 

We consider the following Lagrangian density of a singlet scalar field interacting with 
another scalar field x ^"^^ ^ fermion ip. 

I^=\ {d,<\>f-\ ^J0'-^A 0^+i {d^xf-\ m\x^-M<\>x^-\ h\'<p'+ti,^>^d,^P-m^i,^-f<l)^ij. 

(1) 
When we investigate the time evolution of 0, only the initial condition is fixed, and so 

the time contour in a generating functional starting from the infinite past must run to the 

infinite future without fixing the final condition and come back to the infinite past again. 

The generating functional in the in-in formalism is thus given by 



Z[J,K,r],ri] = Tr 
= Tr 



Tp i exp 

T_ < exp 

X T+ < exp 



P 



dt / d^x{J^(f)+ + K+X+ + V+'^+ + V+i' 



(2) 

where the suffix c represents the closed time contour of integration. X_|_ denotes a field 
component X on the plus-branch (— oo to +oo) and X_ that on the minus-branch (+oo to 
— oo). The symbol Tp represents the time ordering according to the closed time contour, 
namely, T+ the ordinary time ordering, and T_ the anti-time ordering. J, K, and r], fj 
represent the external fields for the scalar and the Dirac fields, respectively. In fact, each 
external field J^{K^,r]^,f]^) and J-{K-,ri-,fj-) is identical, but for technical reasons we 
treat them differently and set J+ = J^{K^ = K_,ri+ = ri^,f]+ = fj) only at the end of 
calculation, p is the initial density matrix. Strictly speaking, we should couple the time 
development of the expectation value of the field with that of the density matrix, which 
is practically impossible. Accordingly we assume that deviation from thermal equilibrium 
is small and use the density matrix corresponding to the finite-temperature state with the 



exception that the low-momentum modes of may have a larger amplitude initially, whose 
fate we are interested in. Then the generating functional is described by the path integral 
as 

Z[J,K,7],r]] = exp(^tW[J,K,r],f^]^ 

Vct> jvxjvi, jvr exp (2^(0, X, ^, ^, J, i^, r], r^ ] ) , (3) 

where the classical action S is given by 

S[(f),x,'ip,'4',J,K,r],f]]= d'^x[C + J{x)(f){x) + K{x)x{x) +r]{x)ip{x) +f]{x)ilj{x)] . (4) 



As with the Euclidean-time formulation, the scalar field is periodic and the Dirac field anti- 
periodic along the imaginary time direction, with 0(t, x) = (j){t—ij3, x), x{t, x) = x(t— i/?, x), 
and ilj{t, x) = —ip{t — i(3, x). Here (3 is the reciprocal of the temperature T. 

The effective action for the scalar field is defined by the connected generating functional 

as 

r[0] = W[J,K,ri,f]] - f d^xJ{x)(l){x) , (5) 

where (f){x) = 6W[J,K,ri,fi\/6J{x). In terms of the components along the plus and the 
minus branches, it reads 

/oo /• 

dt / rf^x [J+(a;)0+(a;) - J-(^)0-(^)] , (6) 

00 J 

with (f)+{x) = 6W[J+, J_, ■ • ySJ+ix) and 0_(a;) = -6W[J+, J_, • ■ •]/5J_(a;). 

We give the finite-temperature propagator before the perturbative expansion. For the 

closed path, the scalar propagator of x has four components consisting of x±{x) ^^^ X±{x')- 



\JT -y I tjy tXj j 



G++ix-x') G+-ix-x') 
G-+{x-x') G--{x-x') 

Tr[T+x+(a;)x+(a;')p] Tt[x-{x')x+{x)p] 
Tr[x_(a;)x+(a;')p] Tr[T^x^{x)x^ix')p] 

G^{x-x') G+{x-x') 
G-{x-x') G^{x-x') 



where 

^ 



G^{k) = — "— + 27fnB{uJk)S{k'^ -m^ 

^ k^ — m^ + te 

Gl{k) = ^, _ ~; _ ^^ + 2vrM^.) ^(fc^ - m;) 

A 



G^{k) = 27r[e{-ko)+nB{uJk)]6{e-ml), 

G~ik) = 2n[e{ko)+nB{uJk)]6{k^-ml), (8) 

with UBiuJk) = {e^"^" - 1)"\ ^fc = Jk^ + ^2 , and e(A;o) = ^(A;o) - e{-ko) ^ . Similar 



formulae apply for (p field as well. 

The propagator for a Dirac fermion is given by 



OquyOC X j 



where 



S:^'^{x-x') 5*^ (x - x' 
5^ + (x-a;') S;p-{x-x') 

Tr [ T+ip+ (a;)^+ (a;')P ] Tr [ —■?/'_ {x')ip+ {x)p ] 
Tr [ ■?/'_ (3^)"^+ {x')p ] Tr [ T_-?/'_ (x)-?/^- (^^Op ] 

5*^(0; -x') 5*^(3; - x') 
S^{x-x') ^J(x-x') 



S'f(fc) = ■ 2nnFiEk){/i: + m^)6{P -ml) , 

^ /k — m^ —It ^ 

S^ik) = 2n[e{-ko)-nF{Ek)]{/c + m^)6{k^-ml), 

S^{k) = 27T[e{ko)-nF{Ek)]{/c + m^)6{e-ml), (10) 



with npiEk) = (e^^^ + l)-\ Ek = J k" + m^ 



■V 



B. Perturbative expansion of the finite-temperature eff"ective action 

The perturbative loop expansion for the effective action F can be obtained by transforming 
~* 0d + C where 0c/ is the field configuration which extremizes the classical action S'[(/), J] 



and ( is small perturbation around (j)ci. Up to two loop order and 0(A^, h'^, p), T is made 
up of the graphs as those depicted in Fig. [T] etc. Summing up these graphs, the effective 
action F becomes 
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^(50+)' - kd^-f - \^l (0+(^) - 0-(^)) - ^ (0t(^) - 0-(:^)) 



4! 






(11) 



where each of i^j [</>+, </>_] corresponds to each graph in Fig. [T]and is given as follows. 



^ 



a 



Li 



e 



Li 



XX 



L2> 







Ls 




L6 




Li 




Ls 



FIG. 1: Feynman diagrams corresponding to each term of the effective action. Solid line denotes 
(J), while broken line and double line represent ip and X; respectively. 



Li 

L2 



-0+(x)G+-(x - x'fcp_{x') + <P-{x)Gj 



'_{x)G-^ + {x-x')\ 
{x-x'f(t).{x')\ , 



[X 



(12) 



(13) 



^3 = ^ / d'xd'x' [<Pl{x)Gl+{x - x')^l{x') - <pl{x)G-+{x - x'fcPlix') 



-^l{x)Gl-{x - x')m{x') + <pl{x)G-f{x - x'f<Pl{x')\ 



(14) 



^f 



d'^xd^x' U+(x)5++(x - x')Sl^{x' - x)(t)+{x') 



-0„(a;)5'7+(x - x')S^'{x' - x)(j)+{x') - (f)+{x)S^'{x - x')5 . +(x' - x)(j)-{x') 



+0_(x)5^~(a; - x')S;j;~{x' - x)0_(x')] , 



' 4 



d'xG++{0)[<Pl{x)-<l)'_{x)] 



(15) 

(16) 
(17) 



L7 = ^ f d^xd^x' [(j)4x)G^+{x ~ x')G++{x - x')^(j)+{x') 
— 0_(x)G. ^(x — x')G'^^(x — x'Y(p+{x') 



-0+(x)G+~(x - x')G^"(x - x')V-(x') 






+0_(x)G^-(a; - x')G--{x - a;')V-(a:')] , 

t/^xrfV [0^(x)Gj+(x - x')V+(a;') - 0^_(x)G^+(x - x')V+(a;') 
-0^(a;)G'+-(x - x')'0^(x') + (j)l{x)G--{x - x')^(j)l{x')] . 



It is convenient to introduce new variables 



1 



(t>c = ^{(f)+ + (1)^) and (j)A = (/)+- (j)- 

to rewrite the effective action in terms of tliese variables. As will be seen in 
response field and (pc is the physical field. We find 



(18) 
(19) 

(20) 

, 0A is a 



r[0c,0A] 



d^x 



A 



-(Pa{x)0(J)^{x) -mj0A(a;)0c(a;) - tt (0i(x)0c(a;) +40A(a;)0c(x)) 



4! 






(21) 



with 



Li 



2 d^a;G++(O)0A(x)0,(x), 
A^ /" 



iX^ 
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+— / d^xd^x'Re [G++(x - x')'] 0A(x)0A(a;'), 



(22) 



(23) 



L^ = -^ j d''xd^x'\m[G++{x-x'f] {M^)<Pc{x)<Pl{x') + \M^)<Pli.^')<Pc{x))e{t,-t2) 

+^ I d^xd^x'Re [G++{x - x'f] 0A(a;)0A(a;')0c(a;)0c(a;'), (24) 
L4 = 2/2 I d'^xd^x'lm [S++{x - x')S^+{x' - x)] <pA{x)<pc{x')6{ti - ts) 

_!/_ I d^xd^x'Re [S++{x - x')S^+{x' - x)] 0a (x) 0a (a:'), (25) 
L5 = -4A<2 f d'^xd'^x'Im [G^^{x - x'Y] (f)A{x)(Pc{x')e{ti - h) 

+iM f d^xd'^x'Re [G^^{x - x'f] 0A(a;)0A(a;'), (26) 

Le = -^ J d''xG++iO)Mx)Mx), (27) 
L7 = -h^ f d^xd'^x'lm [G'++(x - x')G^+ix - x'f] 0A(a;)0c(x')^(ti - ts) 

+ -^ / d^xdVRe [G++(x - x')G++{x - x'f] Mx)Mx'), (28) 
Ls = -y /"d'xciVlm [6-++ (a; - x')'] [0A(a:)0e(x)02(a;') + i0A(a;)0i(x')0e(a:)] 9ih - ^2) 

+ ^— / rf^xdVRe [G'++(a; - x'f] 0A(a;)0A(a;')0c(a;)0c(a;'). (29) 

Among these terms, Li and Lg are corrections to the mass term of 0, while other terms have 
both real and imaginary parts. As a result we find 

r[0e,0A] = fd^x |-0A(a:)[n + M^](j),{x) - ^ [40a(x)03(x) + 0e(x)0i(a;)] 

d'^xd'^x' [A2{x — x') + Ai{x — x) + ^5(3; — a;') + A-j{x — x') ] ^(ti — ^2) 0a (a;) 0c (a;') 

d^xd^x'\A^{x - x') + ^8(3; - x')]e{ti - ta) 

X [0A(a:)0c(a:)0^(a;') + i0A(a;)0i(a;')0c(a;)] 
+- / d^xd^x' { [B2{x - x') + 54(0; - x') + B5{x- x') + Bj{x - x')] <Pa{,x)<Pa{,x') 

+ [B^{x - x') + B^{x - x') ] 0A(a;)0A(a:')0c(x)0c(a;') } , (30) 

where 

M' = ml + \N{m,) + 2h'N{m^\ N{m.) ^ /^^^^^, ^. = ^"^7^,(31) 

A2{x-x') = ^lm[G^{x-x'f] , (32) 



Asix- 


- x') 


Ai{x- 


- x') 


Anix- 


- x') 


A'j{x- 


- x') 


A^{x- 


- x') 


B2{x- 


- x') 


Bsix- 


- x') 


B,{x- 


- x') 


B,{x- 


- x') 


B^{x- 


- x') 


Bsix- 


- x) 



Um[Gl{x-x'f] , 
-2/2lm {tr [S^{x - x')S^{x' - x)] } 
AMHm[Gl{x - x'f] , 
hHm[Gl{x - x')G^{x - x'f] , 

^lu.[Gl{x - x'f] , 

Re[G^{x-x')'] , 

Re[G^ix-xr] , 
-fRe{tT[S^{x-x')S^{x'-x)]} , 
2M^Re[G^{x-x'y] , 



A2 
6 

2 



Y 



Re [G^ix - x')G^ix - x'f] , 
Re[G^{x-x'y] . 



(33) 
(34) 
(35) 
(36) 
(37) 

(38) 

(39) 
(40) 
(41) 
(42) 

(43) 



Apparently, Aj and Bj are related with the real and the imaginary parts of Lj, respectively. 
The above expressions for Aj{x — x') are valid only for ti — t2 > 0. We find 



Aj{x' - X,ti-t2) = Aj{x -x',ti-t2), 
Aj{x — x',t2 — ti) = —Aj{x — x',ti—t2) 



(44) 
(45) 



for t2 — ti < 0, although only those with ti — 12 > appear in the final expressions. We also 
find 



Bj{x' - x,ti -t2) = Bj{x - x',ti -t2), 
Bj(x — x',t2 — ti) = Bj(x — x',ti — t2). 



(46) 
(47) 



The imaginary parts of the effective action represent dissipative effects and we can obtain 
real effective action by introducing auxiliary random Gaussian fields, C,a{x) and C,m{x), as 
follows. 

exp(ir[0c,0A]) = / T^U / 'I^^mPa[^a]Pm[^m]exp{iTcs[(j)c,(l)A,^a,^m]} , (48) 



where 



rcff[0c,0A,Ca,Cm] = Rer[0c,0A] + / d^x[^a{x) 



^a{x) + ^m{x)<l)c{x)<f)A{x) ] ■ 



(49) 
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Here Pa[Ca] and Pm[^m] are the probability distribution functional defined by 

1 



Pam 



ATaexp 



d Xd x'^a{x)D^ (X - X')^a{x') 



Da{x - x') = B2{x - x') + B4{x - x') + B^ix - x') + Br{x - x') 



Dm{x — x') = B^{x — x') + B^{x — x' 



i I d'xd^x'U{x)D;^\x - x')U{x') 



(50) 
(51) 
(52) 
(53) 



respectively. Thus the dispersions of C,a{x) and ^m{x) are given by Bj. In the above expres- 
sions Ma and Mm are normalization factors, while the inverse D^^{x — y) is defined by the 
relation 

' d^yDr\x - y)D,{y - z) = 6{x - z). (54) 



C. Equation of motion 

Applying the variational principle to Feff, we obtain the equation of motion for (^^ con 
taining no imaginary quantity. 



0. 



(55) 



0A=O 



From (Uni), it reads 



( D + M^ ) 0,(x) + \^1{X) + f dt' I d^x'Caix - X')0c(x') 



with 



+0e(x) / dt' / d'x'Cmix - X')<f)iix') = U^) + M^)Uix) , (56) 
J —oo J 

Ca{x - x') = A2{x - x') + A^ix - x') + A5{x - x') + Aj{x - x'), (57) 

Cm{x-x') = A3{x- x') + Asix- x'). (58) 

We shall call these two functions memory kernels because the last two terms in the left-hand- 
side of fl56|) are nonlocal in time. They will reduce to the dissipation terms and perturbative 
corrections to the classical equation of motion which would become a part of the derivative 
of the effective potential, V^g(0), if we restricted (j){x') to be a constant in space and time. 
In this equation of motion ^a{x) should be regarded as an additive random Gaussian noise 
with the dispersion, 

{U^)U^')) = Da{x-x'), (59) 
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and Cm{x) is a multiplicative random Gaussian noise acting on (pdx) with the dispersion, 

{^m{x)^rn{x')) = Dm{x -x') . (60) 



III. ANALYSIS IN THE LINEAR REGIME 



A. Equation of motion in the Fourier space 



Here we concentrate on the case only linear terms of (pc are important and multiplicative 
noise is negligible in the equation of motion (j56p . Then the equation of motion reads 



( D + M^ ) (P^x) + dt' d^x'Caix - x')0c(a;') = ^a{x) . (61) 

J —oo J 

Hereafter we omit the suffix c. 

In this regime it is convenient to rewrite the above equation in the wavenumber space. 
Defining the spatial Fourier transform as 

0,(t) = j d'x<j){x, 1)6-'"-, Ut) = J d'xUx, t)e-"'-, (62) 

Ca,it-t') = fd^xCaix,t~t')e-"'-, (63) 

we find 

Mt) + {k' + M^) 0,(t) + f dt'Ca,{t - t')0.(t') = Ut) , (64) 

where Cak{t — t') is a real function thanks to (j44j) . Here the noise term in the Fourier space, 
^fc(t), is a random Gaussian variable with the dispersion, 

{Ut)C'it')) = j d'xDa{x,t- t')e-"'-{27rf6ik - k') = Da,{t - t'){2Tif5{k - k'). (65) 

Thus each Fourier mode is completely decoupled from each other in the linear regime even 
in the presence of the noise term, as it should be. 

Equation (j64p can be solved in terms of the Fourier transform with respect to t, 

0(^) = j dtct>{t)e'^\ U^) = j dtUt)e'''\ (66) 

Ca,{uj) = j dtCa,{t)e'''\ Da,{uj) = j dtDa,{t)e''^'. (67) 

Here note that Cak{uj) is pure imaginary due to (ji3|) . Using the formula 

/ dre'^"^-^'^^ = iF + 7cS{iu-uj'), (68) 

Jo uj -u' 
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we find 




i-u 


' + k^ + M^)M^) + 


Defining 


real quantities 




m2 = M' + e 



^F^-^tCa,{uj')M^) + ]-Ca,{uj)M^) = U^). (69) 
Itt io — io 2 



^P^_,C'„,(c.O, t,{uj)^z^^^^, (70) 

which respectively constitute real and imaginary parts of the self energy of 0, we obtain 

7 27r ,,2 _ ^2 ^ ,^ffc(^) 7-00 J 27r ,;2 _ j\^2 ^ ,^ffc(o;) ■ 

If Tk{^) satisfies < r^.(u;) ^ M^ and cu-dependent part of M^ is negligibly small, which 
turn out to be the case in the specific examples discussed later, (f7T|) has poles at u 
±Mfc — irk{Mk)/2 and it can be solved as 



rvj 



Mt) = 4r f dt'e-'^^'^^^'^^^'''%mMk{t -t'Mt'). (72) 

Mk y_oo 

Adding two independent homogeneous modes, a general solution with an arbitrary initial 
condition (f)k{ti) and (f)k{ti) at some initial time t = tj is given by 

Mt) = Ut;)e-MM.)it-u) ^^^M,{t - U) + ^e-^f^(^^^)(*-*") sinM,(t - U) 

+ Yr f dt'e-'^^-^^'-^^'-'^smM,{t-t')Ut'). (73) 

^vik Jti 

by virtue of the assumption Tk{Mk) -C M^. 

Then using (jU3j) and (pTjl . the expectation value of the absolute square amplitude at late 

time t > ti + f^^(Mfc) reads 



mm- '''"'^^'^ 



1 H — — sm 2Mut 

Mk 



(27r)3<5(0). (74) 



2MiTk{M, 

The second term in the bracket vanishes of course if we take time average over an oscillation 
period as well. Equations (fTSj) and (fTIjl indicate that each mode does not decay completely 
but its square amplitude approaches an equilibrium value determined by the ratio of the 
power spectrum of the noise to Tk{Mk) with the time scale r^"'^(Mfc) = —2iMk/Cak{Mk). 

In order to evaluate these quantities we must calculate Fourier transform of the memory 
kernel Ca and the noise correlation Da explicitly using the expressions given in the previous 
section. Below we study the effects of interactions with fermions and bosons separately in 
turn, because they have different behaviors due to the different statistical properties. The 
striking difference between fermionic noises and bosonic noises have been pointed out in |l7|. 

13 



B. Interaction with fermions 



First we study the case the scalar field interacts only with fermions. In this case Ca and 
Da are governed by L4, namely (j3^ and (jlOl)- Because both A^i^x — x') and B^^^x — x') 
are parity-even functions, see (j44p and (|46|) . the spatial Fourier transform is identical to the 
Fourier cosine transform, so we find 

Cau{t - t') = A,,{t - t') = -2f J d^xe-^'-lm {tr [S^ix,t- t')S^{-x,t' - t)] } 

= -2flm {tr [5j(p, t - t')S^{p + k,t'- t)] } . (75) 

Here we have used the fermion propagator expressed by the real time t and spatial wavenum- 
ber p, 



S^{p,t) 



d^xS^{x,t)e~'''-^ 
Epjo^P2±m^^p ^-iE,t _ -Eplo -Pl + m^ _ ^F.iE.t 



2E„ 



2E„ 



9{t) (76) 
di-t), 



with Ep = Jp^ + m^ and n^ = np{Ep). Using 



tr[5j(p,r)5,J(p + fc,-r) 

- i-EpE.+p -p.k-p' + mJ)«+pe^(^-+^^+'=)- 

- i-EpE,^p -p.k-p' + ml){l- nf )(1 - nf^^)e-(^^+^^+'=)^ 
HEpE.^p -pk-p' + mj)<(l - nf+^)e^(^-^^+'^)^] , (77) 

with r = t — t' > 0, we find that the Fourier transform of the memory kernel is given by 

d^p 1 



Cak. (uj 



drAAjT)e'- 



-2nxt 



(27r)3 EpEk+p 
{{EpEk+p -p^ + ml) [{1 - n^)n^+p - n^{l - n^+p)] S{uj + Ek+p - Ep) 

+ {EpEk+p -p^ + ml) [nj(l - n^^p) - (1 - n^)n^+p] S{uj - Ek+p + Ep) 

+ {EpEk+p + p'- ml) [n^n^^p - (1 - <)(! - n^^p)] 6{uj + E^+p + Ep) 

+ {EpEk+p + p'- ml) [{1 - 0(1 - n^^p) - nf nf J 6{uj - E^+p - Ep)} , (78) 
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where -Efc+p = \ {k + pY + m^ and ^f+p = np^Ek+p). The first term in each square bracket 
can be interpreted as decay or absorption of 4>k{^), which is denoted by Re, while the 
second term corresponds to inverse decay or creation of (pk{^) denoted by Re, because n^ 
corresponds to the number density of an initial state and 1 — nf^ to the Pauli-blo eking factor 
of a final state. The above expression (f7H|) is closely related with the discontinuity of the 
self energy of (p at finite temperature which was obtained by Weldon J2^ using a different 
procedure. In his approach one had to add and subtract appropriate combinations of rip and 
"^k+p to obtain the above form in which physical interpretation of absorption and creation 
of is manifest, while in our scheme the above result is obtained straightforwardly from the 
structure of the fermion propagator (fTUI) . 

Due to the delta function the ratio of creation and destruction rates satisfies the detailed- 
balance relation, 

^ = e-'-, (79) 

for all combinations. For example, in the first square bracket of the right-hand-side of (|7H|l . 
we find 



(80) 



R£ ^ np (1 rifc^p) ^ ^^^^^^_^^^ ^ ^^^ 
Rd (1 - <)nf+, 

under the condition u = Ep — E^+p coming from the delta function 6{uj + -E^+p — Ep). 
The dispersion of the stochastic noise in Fourier space, on the other hand, reads 

Dfuj) = I dTB^JT)e'^^ = nf f -^-^ — 
^ ^ y ^ ' J {2TifEpEk+p 

X {{EpEk+p -p^ + ml) [{1 - nl)nl^p + <(1 - n^^p)] 6{uj + E^+p - Ep) 
+ {EpEk+p -p' + ml) [nf (1 - n^^p) + (1 - Onf J 5(^ - E^^p + Ep) 

+ {EpEk+p +p^- ml) [nX+p + (1 - <)(! " ^k+p)] ^(^ + E^+p + Ep) 

+ {EpEk+p +p'- ml) [{1 - 0(1 - n^^p) + <nf J 6{uj - Ej^+p - Ep)} , (81) 

In this dispersion, both destruction Rd and creation Re contribute additive manner. From 

dZOl), dZHl) and dHU), we find 

£^ = -2,.aiM=.^£±^=.f^^2T. (82) 

r*(i^) C„(w) fie --Re e»^-l 

where the last approximate equality holds for the soft modes with uj <^T. This is nothing 

but the fluctuation-dissipation relation derived purely from quasi-nonequilibrium quantum 

field theory at finite temperature. 
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Note that the fluctuation-dissipation relation has also been obtained by Gleiser and 
Ramos |l5[ in the context of nonequilibrium fleld theory at flnite temperature. However, 
because they assumed the scalar fleld evolves adiabatically, they had to invoke higher loop 
effects to obtain a nonvanishing dissipation coefficient. As a result their noise term and dis- 
sipation term appear at different order of perturbation. This problem of adiabatic treatment 



has been pointed out by Gleiner and Miiller 18| who adopted a harmonic approximation 



instead in order to extract a term proportional to 0, which represents dissipation in the 
equation of motion and obtained the correct result. In the present analysis we have made 
no assumption about the adiabaticity of the evolution of the scalar field but worked in the 
Fourier space assuming that the quartic term dominates its potential. Then we can see that 
both dissipation and noise terms appear at the same order of perturbation. 
The time average of ()74|) over an oscillation period reads 

l{\Mt)\') = lM',{\Mt)\') = ln2n)'6{0). (83) 

This equation shows the classical equipartition law is satisfied for low-momentum modes 
that kinetic energy per degree of freedom is equal to T/2. This property can be seen more 
manifestly if we adopt a box normalization with a finite side L and periodic boundary 
condition. Then (f){x, t) is expanded as 

0(^,t) = 5^0„(t)e^^"-, (84) 

n 

with n being a spatial vector consisting of integers. Then (jU3j) is replaced by 

{Ut)CAt')) = Dar.it - t')6^^,, DUt -t') = j ^Daix, t - t')e-^^"-. 
From ()74j) we find average kinetic energy of each soft mode is given by 



(85) 



lilUm = ^-T. (86) 

The above is the results for the Rayleigh- Jeans regime u <^ T, where classical analysis 
applies. We now consider a more general case. Instead of taking the high temperature limit 

u/T — ^ as in the last equality of (jH^ . we rewrite (|H^ as 

D,,iu) Rp + Rc e^^ + l ( , .l\ ,„„. 



tk{u) Rd-Rc e/5--l V 2 
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Then (fTIjl reads 






V 



(88) 



where V = {27t)^6{0) denotes (infinite) spatial volume. Its interpretation is obvious. The 
left-hand-side represents energy density stored in the fe-mode and the right-hand-side shows 
it consists of thermal and vacuum quanta with energy level Mk in the final equilibrium 
state. Thus the interaction with a thermal bath drives each Fourier mode (p^it) to the 
thermal equilibrium value with the same temperature in the time scale rfc(Mfc)~^. 

Next we evaluate the dissipation rate using (fTHj) . Since we are primarily interested in the 
fate of the homogeneous coherent mode, we take k = 0. Then only the last term of ()78p is 
nonvanishing and we find 



r^(T)=fo(Mo 



87r 



Mo 



2m^ 
1^ 



3/2 



l-2nF 



Mn 



r^(o) 



2n, 



Mn 



(89) 

We thus find the dissipation rate at finite temperature is suppressed by the last factor in (|H^ 
due to Pauli blocking with ri7'(0) being the decay rate of a particle at rest into two fermions 
ip and ip at zero temperature. Note that the above dissipation rate vanishes when Mq < 2m^. 
In this case the coherent oscillation is not thermalized through the Yukawa interaction at 
one-loop, because not only the dissipation kernel Cao{uj) but also noise correlation Dao{uj) 
vanishes in this case since both contain delta functions with the same arguments. 

The above arguments are based on the propagator ()76|) where only the zero-temperature 
intrinsic mass m^ is taken into account. If the Yukawa interaction fipipip generates large 
oscillating mass to ip, decay of into two fermions would be possible only during a short 
interval when f\(j)\ < Mo/2 as the scalar field passes through the origin twice in each os- 
cillation period. The dissipation rate of the scalar field in such a situation cannot be dealt 
with the perturbation theory we are using. This issue has been investigated by Dolgov 



and Kirilova 



25[ using a quasiclassical approximation at zero temperature. They find that 



the dissipation rate of (p is not exponentially suppressed but by a factor ~ [MQ/m^ ___V/'^ 



"XjOscy 



where m^^osc is the maximum of x's mass with the oscillating component taken into account. 
Finite-temperature generalization of the analysis in such a regime is not straightforward and 
we restrict our analysis to the perturbative regime f\4>\ ^ Mi here. 

On the other hand, recently Kolb, Notari, and Riotto J26| argue that if the would-be 
decay products of the oscillating infiaton scalar field acquire a thermal mass larger than 
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the inflaton mass in the thermal background, the inflaton cannot decay into these particles, 
and that reheating is suspended for some time based on the observation that the phase 
space would be closed for the mass of the decay product being larger than half the inflaton 
mass. This phenomenon could be observed in our perturbative approach as well, if we used, 
instead of finite-temperature "bare" propagator (f7B|). the "dressed" propagator in which 
finite-temperature higher-order quantum corrections are taken into account. But use of 
such an dressed propagator can easily result in overcounting of diagrams and the detailed 
comparison of two different expansion method is still under way. Here we focus on the 
effects from lowest possible orders and continue to use the bare propagators. In the practical 
applications in §V]and ^Vll we mostly consider the cases the thermal mass of decay products 
remain smaller than the angular frequency of the oscillating field, so both approaches give 
the same results. 

C. Interaction with bosons 

Next we consider the effect of three-body interaction Ai(f)X^ for which Ca and Da are 
determined by L5, namely (j35|) and (|4ip. Again their spatial Fourier transform is identical 
to the Fourier cosine transform due to the parity evenness, and the memory kernel reads 

Ca^ir) = A,,iT) = AM' J d^xe-'''nm[G^{x,TY] 

= ^M' J ^Im [Gj(p, r)Gj(fc - p, r)] , (90) 

for r = t — t' > and Cakir) = —Cak{.—T~) for r < 0. Here G^{p, r) is defined by 



Gj(p,r) = j d\G 



J(a;,r)e--- 



1 



2ujp 



{ [1 + riBioOp)] e-'^-\^\ + nB(^p)e*"-l^l} , u, = ^/p^ + ^2 . (91) 
We obtain 



(27r)3 UpUk-p 
X {[(1 + np)(l + Uk-p) - UpUk-p] 6{uj -Up- Uk-p) 

+ [(1 + np)nk-p - (1 + nk-p)np] 5{uj -Up + ujk-p) 

+ [np{l + Hk-p) - (1 + np)nk-p] 6{uj + Up - ujk^p) 

+ [upUk-p - (1 + np){l + Uk-p)] 5{uj + ujp + uJk-p)} , (92) 
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where Uk-p = w(A; — p)2 + m^, Up = nsi^jp), and rik-p = nB{uJk-p), respectively. The 
first term in each bracket represents destruction {Rd) while the second term corresponds to 
creation (-Re)- Due to the delta function their ratio satisfies the detailed-balance relation, 
Rc/Rd = e~^^ 1 for all combinations. We can also confirm that Tk{oj) = iCak{^)/'2uj is 
positive definite. 

Next we consider the power spectrum of thermal noise given by 



Da.{r) = B,,{r) = 2M' J ^Re [Gj(p, r)Gj(fc -p,r 



(93) 



Its Fourier transform with respect r reads 



Dak{uj) 



n 



-M' 



d^p 



2' ' J (2-kY ujpUJk-p 

X {[(1 + ?2p)(l + Hk-p) + UpTlk-p] 5{UJ -Up- Uk-p) 

+ [(1 + np)nk-p + (1 + nk-p)np] 5{uj -ujp + cUfc-p) 
+ [np(l + Uk-p) + (1 + np)nk-p] S{u + Up - LOk-p) 
+ [upHk^p + (1 + np){l + Hk^p)] 5{uj + ujp + ujk-p)} , (94) 

As in the case of Yukawa interaction (jHT|) . delta functions in ()94j) have been multiplied by 
Rd + Rci which means that both destruction and creation act as a noise to the evolution 
of the scalar field in the same way. 

From ()92j) and (jMj) we find again that 



Dakjuj) 



Rd + Ri 



LU- 



c 



3^"^+ 1 



UJ- 



naniely, 



Rd - Re " e^- 

Dak{Mk) 



2u ( nBioj) 



\\, 



2T, 



(95) 



(96) 



Tk{Mk) 
in the Rayleigh- Jeans limit Mk <^ T . Thus the fiuctuation-dissipation theorem is satisfied 

in this case, too, and the final equilibrium configuration has the same property as in the 

case thermalization proceeds through Yukawa interaction. 

As explained in the previous section, the dissipation rate toward thermal equilibrium 

distribution is given by Tk{Mk). For the coherent zero- mode, in which we are primarily 

interested, one can easily find 



TB{T)=f^{Mo) 






1 



2m, 

Afn 



1/2 



l + 2nB 



Mn 



^b{0) 



l + 2nB 



Mn 



(97) 
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because only the first S function in ()92j) gives nonvanishing contribution when k = 0. Again 
rB(0) is the decay rate of into two x particles through trilinear interaction. Thus the 
dissipation rate is enhanced at finite temperature due to the presence of bosons. In the high 
temperature limit j3Mq ^ 1 (j97|) reads 

1/2 

(98) 



4T M^T 

Ts{T) ^ —Ts{0) = ^ 



Mn ' ' 2nM^ 



1 - ' ^ 



Mo 



Note that these dissipation rates vanish when Mq < 2?7i^. In this case the coherent oscillation 
is not thermalized through one-loop of the three body interaction Ai(j)x'^- For the same 
reason described in the latter part of ^IIIBl the above dissipation rate applies only for 
Mq > Ai(f). In the large field-amplitude regime when this inequality is not satisfied, particle 
creation through broad parametric resonance would be much more efficient |32J|. 

D. Setting-sun diagrams 

Next we study the contribution from the setting-sun diagrams, L2 and L7. Since L7 is 
expected to give larger contribution we first analyze the Fourier transform of the memory 
kernel corresponding to it, which is given by 



Calico) = I dtA-r,iT)e'^^ 



h' J dre^-^ J -^.i^^MPi +P2+P3- fc)Im Mip,, r)G^iP2, r)G^4P3^ r)] 

J {^T^) SUJ1UJ2UJ3 

X { [(1 + ni){l + n2){l + ^3) - nin2n:i] 6{uj - uJi - UJ2 - 1^3) 
+ [nin2n3 - (1 + ni){l + n2)(l + n^)] 6{uj + cui + 6^2 + UJ3) 
+ [ni{l + n2)(l + n3) - (1 + ^i)^2'^3] ^i^ + lui - LU2 - UJ3) 
+ [(1 + ^1)^2^3 - ni{l + n2)(l + ^3)] ^(^ - LU1 + UJ2 + UJ3) 
+ [(1 + ni)n2(l + ris) - ni{l + ^2)^23] 6{uj - coi + UJ2 - u;3) 
+ [ni{l + ^2)723 - (1 + 721)722(1 + 713)] 6{uj + uji- UJ2 + UJ3) 

+ [711^2(1 + r7.3) - (1 + 72i)(l + ^2)^3] 6{UJ + UJ1+UJ2- UJ3) 

+ [(1 + ni){l + 712)713 - 7iir7.2(l + %)] 6{uj -UJ1-UJ2 + UJ3)} . (99) 
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Here we have defined tui = \/pi + rn^, L02 = \ P2 + "^x' ^^ ^ V -Ps ~'~ "^0' ^^"^ ^* ~ '"'^('^i)- 
Once again the first term in each coefficient of delta functions represents destruction {Rd) 
while the second term corresponds to creation {Re)- Due to the delta function their ratio 
satisfies the detailed-balance relation, Re/ Rd = e~^'^, for all combinations as before. We 
can also confirm that Tk{u!) = iCakiy)/'^^ is positive definite. 
Similarly, the Fourier transform of noise correlation reads 



Dak (uj) 



dtB7k(T)e 
2 



iujt 



I "^^"^"^ / Wf^^^'^'^^P' +P2 + P3- fc)Re [G^{p„ r)G^{p„ r)Gl{p,, r 



{2n f 6 {p^+p^ + p.,- k)- 



2 J (27r)9' "' ''^' ■ -^ ■ -^ "'8u;iUJ2UJs 
X { [(1 + ni){l + n2){l + ns) + nin2ns] 6{uj -UJ1-U2- u^) 



+ [nin2n^ + (1 + ^i)(l + '^2)(1 + ^3) 
+ [ni(l + n2){l + ^3) + (1 + ni)n2n'i 
+ [(1 + ni)n2n^ + ?2i(l + n2){l + ^3) 
+ [(1 + ni)n2(l + n^i) + ni(l + n2)n3 
+ [ni(l + ^2)^3 + (1 + ni)n2{l + ^3) 
+ [^1^2(1 + ^3) + (1 + ni){l + ^2)^3 
+ [(1 + ?2i)(l + ^2)^3 + 721^2(1 + ^3) 



5{uo + UJ1+UJ2 + to'3) 
5{uj + uji— UJ2 — (^3) 

5{UJ — UJ1+UJ2+ UJ^) 
5{UJ — UJi+ UJ2— UJ^) 
5{UJ + UJi— UJ2+ UJz) 
5{UJ + UJ1+UJ2— UJ^) 

5{uj — coi — UJ2 + ^3)} • 



(100) 



Each coefficient of delta functions consists of Rn + Rq as before. One can also calculate the 
respective quantities for the other setting-sun diagram L2, which has also been calculated 
in 1^, by the following replacement: 



h' 



AV3, 



Gj(Pi, r) -. G^Jp,, r), ^^(p^, r) ^ G^Jp,, r) 



(101) 



^i — ^ a/Pi + "^i ^2 



pI + iv?., with Hi 



nsiuJi). 



In both cases we find the same structure again for Dak{'jj) and f a;('^) = iCakiyj)/'^^, that 



IS, 



bak{uj) Rd + Rc e^^ + i / . . r 

— = UJ = tu— 5 = 2uj nniuj) H — I , 

TJu) Rd-Rc e^^-l V 2 " 



(102) 



and the fluctuation-dissipation relation is satisfied. 
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Since the setting-sun diagrams involve three particles in the intermediate states their 
contribution for the zero mode is nonvanishing even when Mo/3 is smaller than the mass of 
the interchanged particles. Hence this two-loop effect could be important when Mq is smaller 
than 2m^ or 2r7i^ so that one-loop effects discussed in the previous sections are inoperative. 

The analytic evaluation of the dissipation rate with these diagram is cumbersome for 
general cases, so we report the result only for several limiting cases. First in the high 
temperature limit with T ^ Mq ^ m^, m^, it reads 

for A(/)'^/4! interaction, and 

for /i^ 0^x^/4 interaction, where we have used a formula 

"^ , Inx TT^ , , 

and neglected m^ and m^ in the intermediate state. In this case, the first, the third, the 
fifth, and the last delta functions of ()99|) and ()100j) give nonvanishing contribution. 

Since h^ is likely to be of order of A (^ A^), we expect (jl(J4p is much larger than (jl(J3|) . 
so we concentrate on the diagram L7 with /i^ 0^x^/4 interaction from now on and take the 
masses of interchanged particles into account. Then we find 

for T ^ Mq = m^ ^ m^. In this case the third, the fifth, and the last delta functions of 
(I99|) and (jlUUj) give nonvanishing contribution. Even in the case the mass of interchanged 
particle m^ is much heavier than Mq one can easily see that 6{Mo — uji + uj2 — UJ3) and 
S{Mq -\- uji — UJ2 — UJ3) in (j^J^ and ()100|) can give nonvanishing contributions because the 
large masses in ui and ci;2 tend to cancel each other in these delta functions. As a result we 
find that, contrary to the case of Yukawa coupling and three-body bosonic interaction, the 
dissipation rate due to the setting-sun diagram is nonvanishing even if T > m^ ^ Mq > m^, 
and reads 

Thus it is suppressed only by a factor Mq/ui^. This is the case m^ is large and constant. We 
have not manipulated the case m^ has a large oscillating component, but the suppression 
might be even milder then. 
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E. Summary of this section 

Here we summarize the results of our analysis for the case obeys linear equation of 
motion with an additive stochastic noise term (j61|) . Working in the Fourier space we have 
solved the equation and obtained a general solution (f7H|) whose damping rate is proportional 
to the imaginary part of the Fourier transform of the memory kernel, Cak{i^), related to the 
self energy. 

We have then shown that the expectation value of the square amplitude of each wavenum- 
ber mode relaxes to a specific value determined by the ratio of the Fourier transform of the 
dispersion of noise correlation to the dissipation rate. In the Rayleigh- Jeans limit this ratio 
reduces to the temperature and the classical fluctuation-dissipation theorem holds there. In 
more general cases we find that the energy density of each mode is the sum of thermal and 
zero-point vacuum contributions in the final equilibrium state. These results are entirely due 
to the detailed balance relation ()79|) and is independent of the nature of the intermediate 
state in the loop diagram. 

On the other hand, the high-temperature behaviors of the dissipation rate is totally 
different whether it arises from fermionic interaction or bosonic interaction. Although it 
is equal to the decay width of particle at zero temperature in the perturbative regime, 
at finite temperature it is suppressed in the former case due to Pauli blocking ()89|) and 
enhanced in the latter case due to the induced effects (jnZj)- These properties have also been 



obtained using a different technique 2^. In both cases one- loop effects are shut off for 
relaxation of the zero- mode field oscillation when Mq is smaller than 2m^ or 2m^. The 
two-loop diagram would be very important in such a situation. The dissipation rate due to 
two-loop setting-sun diagram Lj is summarized as 



^s{T) 






2567rMo 



for T > Mo > m^, m^ 

for T > Mo = m^ > m^, (108) 

for T > m^ ^ Mo > m^. 
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IV. NONLINEAR REGIME: EFFECTS OF MULTIPLICATIVE NOISE AND DIS- 
SIPATION 

So far we have analyzed the hnear regime when analytic solution can be found for each 
Fourier mode. Next we step forward to include nonlinear interactions and analyze the effects 
of multiplicative noise and dissipation. To do this we have to deal with the full equation of 
motion (j3Bj) which is not soluble analytically. Hence we can at best hope to extract a term 
representing dissipation in the equation of motion and compare it with the noise correlation. 

In the previous case of the linear equation of motion we were able to find dissipation 
rate without explicitly extracting a terrn proportional to 0, which typically represents dis- 
sipation, unlike previous literatures 15|, |l8[ , because we have solved the equation of motion 
analytically and read off the dissipation rate from the solution. Alternatively, however, we 
may also extract a dissipation term with the correct magnitude in the equation of motion 
without knowing a solution. Here we first describe such a procedure for the linear equation 
of motion as a practice to treat multiplicative noise and dissipation. 

A. Alternative derivation of a dissipation term in the linear equation of motion 

Here we return to the linear equation of motion (|61|) and denote the primitive function 
of Ca{x,t) with respect to t by Ea{x,t), namely, 

Ca{x,t) = -^^Eaix^t), (109) 

or 

Ca{x -X',t- t') = —^Ea{x -X\t- t'). (110) 

Then after integration by parts with respect to t', Eq. ()61|) reads 

{a + M^)(l){x)- I £x'Ea{x-x',0)(l){x',t)+ I dt' j d^x'Ea{x-x')^{x')=ia{x), (HI) 

where we have neglected a contribution from infinite past. It is evident from (j70|) and (jl09|) 
that the Fourier transform of the new kernel, Eak{uj), is related to Tk{uj) as 

Ea,{uj) = 2Tkiuj). (112) 
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Let us consider spatially homogeneous field configuration 0(a;, t) = (pit) or fc = mode, for 
which ()111|) reads 

4>{t)+M''(j){t)-Eaoit = 0)(l){t)+ EaoiT)4>{t-T)dT = Ut), (US) 

where Eao{T) and ^oit) are spatial Fourier transform of the respective quantities with k = 0. 
Although the last term of the left-hand-side of ()113p represents dissipation formally, its 
magnitude depends how the scalar field evolves. For example, if we employ an adiabatic 
approximation such as (pit) = (piti) + (f>iti)it — ti) at this stage, 0(t — r) should be treated 
as a constant in the integrand in ()113|) . In this case the dissipation term vanishes, 

/ Eao{T)(P{t - T)dT = Eao{T)dT(j)=-Eao{uJ = 0)(P = to{uJ = 0)^ = 0. (114) 

Jo Jo ^ 



This is the very reason the previous approach in the literatures |12l . Il5| had to invoke higher- 
loop effects or adopt a different method J28| to yield a nonvanishing dissipation rate in the 
adiabatic regime. On the other hand, if we adopt harmonic expansion around t, 

0(t - r) = 0(t) cos Mr - ^ sin Mr, (115) 

as was done by Gleiner and Miiller |l^, the term proportional to (pit) reads 

/•OO /'OO -| 

/ Ea,{T)ip{t-T)dT^ \ Eaoir) COS MrdT^it) = -Eao{M)(P{t) =fo{M)^{t), (116) 

Jo Jo ^ 

which agrees with our result in the preceding section that has been obtained more straight- 
forwardly. Thus we can see that the Fourier transform of the new kernel Ea{x) gives the 
dissipation rate for oscillating fields even if we did not know its exact solution. 

B. fluctuation dissipation theorem for multiplicative noise and dissipation 

We now apply the above observation for the full nonlinear evolution equation (jSUj) without 
solving it. Since we have fully clarified the roles of additive noise and the corresponding 
kernel Ca{x — x') in the equation of motion in ^1111 we omit these terms and consider the 
following equation. 

( D + M^ ) 0(x) + ^0^(x) + (p{x) I dt' I d^x'Cmix - x')0^(x') = (p{x)U{x) ■ (117) 
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In order to follow the same procedure as in ^IVAI we define the primitive function of Cm{x, t) , 
Em{x, t), in terms of 

CUx,t) = -E^{x,t), (118) 

and perform integration by parts with respect to t', to yield 

( D + M^ ) 0(x) + ^(f>Hx) - I £x'(P{x)Era{x - x\ Q)(P\x\ t) 

+ 1 dt' I d^x'2(l){x)Em{x-x')(f){x')(j){x') = (f){x)U{x)- (119) 



Here the last term in the left-hand-side includes effects of dissipation. So in order to see if 
the fluctuation-dissipation relation also holds for the case of multiplicative noise, we perform 
the Fourier transform of S{x,x') = 2(f){x)Em{x — x')(f){x') as 



S^,,,iuj,uj') = I d-^xdt / d-^x'dt'£{x,x')e-"-^+''^'e- 
d^K dVL^~ 



J 72^^'^^'-''^'' ~ ^)^rnKmk'+Kiuj' + fi), (120) 



where 



EmA^) = f d'xdtEUx)e-'''-^'''\ (121) 

Equation (J12(J|) should be compared with the Fourier transform of the two-point correlation 
of the multiphcative noise, V{x,x') = 4>{x){C,m{x)C,m{x'))(j){x') = (f){x)Dm{x~x')(j){x'), which 
reads 



Vky{uj,uj') = d^xdt d^x'dt'V{x,x')e~ 

r d^K dVt ~ - - 

= —-^—(t)u-K{oo-n)D„,An)(i),.+Auj' + n). (122) 



Here DmK{^) is defined in the same way as (J12H) . 

The multiplicative noise and dissipation under consideration are generated by two graphs, 
L3 and Lg, in the effective action. Since the relevant kernels Ai{x — x') and Bi{x — x') 
{i = 3, 8), namely (jH^ . (jHTj) . (jH^ . and PHjl . have the same structure as in the case of three 
body bosonic interaction, L5 or (J35|l and (J4T| . which has been discussed in §111 (^| we can 
easily obtain the Fourier transform of the kernels. 

Let us first consider the contribution of L3 to CmK{^) = —i^EmKi^) and DmK{^)- From 
(P|) and dni we find 
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and 



A2 [ d^q 1 

= —ill — / - — -^7 

8 J {2ny^ UgUK-g 

X {[(1 + ng){l + riK-q) - nqTlK^q] 5{Vt -UJg- UJK-q) 
+ [(1 + nq)nK~q - (1 + nK-q)nq\ 6{Vl -UJq + LdK-q) 
+ [nq{l + riK-q) - (1 + nq)nK-q] S{Q + UJq- LdK-q) 
+ [nqUK^q - (1 + nq){l + TiK-q)] 5{Vt + UJq+ UJK-q)] , (123) 

X {[(1 + ng)(l + riK-q) + UqUK-q] 5{VL - UJq - tOK-q) 
+ [(1 + nq)nK-q + (1 + nK-q)nq] 6{Q -UJq + UJK-q) 

+ [ng{l + riK-q) + (1 + ng)nK-q] 5{n + ujg - ujK-q) 

+ [UqUK-q + (1 + ^g)(l + TlK-q)] 5{Vt + UJq + UJK-q)} , (124) 



where Uq = J q^ + m^, t^;^„g = ^(K - q)"^ + m|, n, = UB^ujg), and n^_g = nsiujK-q)- We 
can read off destruction terms Rd and creation terms Re of as in ^III CI with Re/ Rd = 
g-/3Q_ Hence we obtain 

DmKJn) ^ 2Q D^^{Q) ^ 2f]^!!!±i (125) 

This is twice the corresponding results for additive noises, (jHH), (P3j) . and ()102p . But this 
discrepancy is compensated by an additional factor 2 in p20|) . so we can see that the 
same relation holds between the noise dispersion and the actual dissipation rate for the 
multiplicative case as in the case of additive noises. Hence the generalized fluctuation- 
dissipation relation is satisfied in this case, too, to establish thermal equilibrium in the final 
state. 

Note that contribution of Lg can also be manipulated by the replacement A^ — > h^, 



UJq — > Jq^ + ml, and ujK-q — ^ \ {K - q)^ + m?,. 
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C. Dissipation rate of zero- mode oscillation 

So far we have shown that multiphcative noise and dissipation also satisfy the desired 
fluctuation-dissipation relation generically, but we cannot obtain the magnitude of the dissi- 
pation rate without knowing the Fourier transform of the scalar field itself. Here we consider 
a specific field evolution and calculate the dissipation rate. To do this we consider the case 
only zero-mode oscillation with a fixed angular frequency M is present, namely, we adopt 
the harmonic expansion as in ()115|) . 



[x, t-T) = (t){t-T)= (f){t) COS Mr - -^ sin Mr, 



(126) 



which is the case we are most interested in. For this approximate solution to be valid we 
assume that A(/)^/3! is smaller than the mass term M^0 in the equation of motion ()119j] . 
Then the last term of the left-hand-side of Eq. ()119|) . which represents the dissipative effects, 
reads 



dt' I d-^x'2(j){x)Emix - x')(j){x') = / dr j ^2E„o(cu)0(t)0(t - r)(j){t - r) 



■OO 

dT I ^2E„o(^)e-*"Xt) 



2tt 



27r 

where 

Ai 



(p{t)^{t) cos2Mt - - I ^-^ - McP^t) I sin2Mr 



M 



- M<p\t) 4>{t) 



(127) 



dr [ ^2E^Juj)e-''^^ cos2Mt = - 



X" 



647rM 



A2 



dr 



27r 



\M ) 

duo ~ 



E^,{2M) + Emo{-2M) 



:i28) 



.1/2 



[l + 2nB(M)] 



27r' 



-EmJuj)e 



sin 2Mt. 



Multiplying the effective equation of motion. 



4>{t) + MV(i) + Ai</>'(t)</)(t) + A2 



by 4>(t), we find 

dt dt\2^^' 2 ^^' 



M 



-M</)2(t)U(t) = 0, 



(129) 



(130) 



-K,<P\t)^\t) - A, 



M 



-M02(t)U(t)0(t). 



;i3i) 
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So far (f){t) is a value at an arbitrary time t around which the harmonic expansion 
is performed. The right-hand-side of the above equation severely depends on the phase of 
the scalar field at the time t. Hence we take an average over the phase of the oscillation to 
obtain its typical magnitude. As a result we find the second term vanishes and (J13H) reads 

^ = -iAi0^(t)p,(t), (132) 

where 0^(t) should now be interpreted as a mean square amplitude around the time t rather 
than its instantaneous value then. 

Thus the dissipation rate is given by 



A^0^(t) 






1/2 



[l + 2nB{M)]. (133) 



1287rM 

Similarly, the dissipation rate associated with h?(j)^x^ 1^ interaction represented by the 
graph Lg reads 



Although the interaction h'^clP'x^ /^ represents creation of a pair of x from pair annihilation 
of formally, the coherent nature of field oscillation makes it possible to interpret the above 
dissipation rate just as a decay of cj) particle with oscillating frequency 2M through three- 
body bosonic interaction M.(t)X^ with the coupling strength M. = h'^(f)/2. They are valid for 
A02 < M2 and h\(j)\ < M. 

V. APPLICATION TO THE LATE REHEATING PHASE OF THE INFLATION- 
ARY UNIVERSE 

So far we have studied relaxation of an oscillating scalar field through various interaction 
channels and obtained the dissipation rate or the relaxation time scale to thermal equilibrium 
for each case. We now apply our results to two cosmological problems, one the reheating after 



inflation j^ |5| and the other evaporation of oscillating quasifiat direction in supersymmetric 
theory in relation with Affleck-Dine baryogenesis [6] . In this section we consider the former 
problem and the latter will be discussed in the next section. 
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A. Brief review of previous results 



First we list several usefu 
extensively in the literatures 



formulae of reheating after inflation which has been studied 



m, 



30, 



31 



Slow-roll inflation is terminated as it is followed 
by coherent scalar field oscillation, whose energy density is released to that of radiation 
subsequently. Two mechanisms are known to reheat the universe. One, which can be very 
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As its name tells, however, 



efficient, is parametric resonance dubbed preheating f 3lL 
it is effective only in the early stage of reheating when the inflaton scalar field is oscillating 
with a sufficiently large amplitude and only when it is coupled to other scalar fields. The 
other is perturbative decay of the infiaton field which terminates reheating process. One 
loop calculation shows that the dissipation rate is equal to the decay rate of the infiaton 




For a constant value of the decay rate F,^, the energy density of the oscillating infiaton, 
P(f,{t), and that of radiation, pr{t), satisfy the following transfer equations. 

dp^(t) 



dt 
dprjt) 

dt 



-{3H + T^)p^{t), (135) 

-AHpr{t) + T^p^{t), (136) 



which are valid when parametric resonance is unimportant. We are also assuming that the 
scalar field oscillation is driven by its mass term and higher order interactions are negligible, 
namely, 

p,(t) = l02(t) + iMV(t). (137) 

The solution of p35|l and p3(i|l are then given by 

P*(«) = P.(«.)(^)'e-^*"-"', (138) 

"'W = "'('■) (^)"^'*/.'(^)"^*<^**- •''" 

Here tj is the time when parametric resonance becomes no longer effective or the epoch 
when the inflaton starts coherent oscillation after inflation, whichever comes later. In the 
latter case we take pr{ti) = of course. In the above system the scalar fleld decays around 
t ^ FT^ and reheating is completed. For deflniteness we deflne the reheating epoch by the 
time when the Hubble parameter H becomes equal to F^, so that the reheat temperature. 
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Tr, reads 



1/4 



/ on \ v^ 
Tn=i^) ^/TW^ = 0.46(f,-i/^ V'TVWg, (140) 

where M^ = 2.4 x lO^^GeV is the reduced Planck mass, and (?* is the effective number of 
the relativistic degrees of freedom with g^, = g^/200. 

Note, however, that this is not the maximum temperature after inflation but that when 
entropy production from the inflaton is practically terminated. Even when preheating is 
inoperative, the maximum temperature can be much higher than Tr as we can write (jl39p 
as 

Prit) = h^t (^^^ p^{U) = ^r^iJM^, (141) 

for tj <^ t <^ Ft^ with Pr(ti) = 0. That is, if the decay product of inflaton is rapidly 
thermalized, the cosmic temperature in the field oscillation regime without preheating is 
given by 

T ^ V^HMl . (142) 



Note that this expression is valid well until the reheating time H = T^ when fll42|) agrees 
with ()140j) with an error of 26%. From ()14nj) and ()142j) we obtain a formula 

T = 0.54^-,-^/® (T^iJMG)^/^ ~ {TlHMcf/^, (143) 

which will be useful later. 

The above is the case with a constant F^. We now consider the cases dissipation rate of 
the inflaton is given by our new results with possible temperature dependence. As in (^ we 
take the interaction Lagrangian as 

- /:int = \ M202 + ^A 0^ + ^ mlx^ + M<Px^ + \ h\^^ + m^V'V' + /0^^ , (144) 

where we assume m^ and m^ are much smaller than the inflaton mass, M, and neglect 
them in the subsequent discussion. As before, M^ includes both intrinsic mass m?, and 
high-temperature corrections of order of ~ h?T'^ and/or ~ /^T^. These thermal masses 
are present if the oscillating masses of x ^i-nd ip are smaller than the temperature, namely, 
/i|0| < T and /|0| < T, respectively. We therefore find K^T'^cj)'^ < T^ and /^T^^ < Tl 
Since the energy density of oscillating inflaton remains larger than that of radiation up to 
the reheating time, these inequalities mean thermal masses are smaller than m^, so M = m^ 
in this regime. 
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It should be understood that the above form of the interaction Lagrangian is a result of 
expansion around the potential minimum which we have set to = after an appropriate 
shift, if necessary. Hence the values of the parameters may not be fixed by the amplitude 
and spectrum of density fluctuations straightforwardly. In particular, if inflation occurred 
more than once, the parameters of the last inflation may entirely be free from large-scale 
observations. If, on the other hand, it describes the original potential as it is and if chaotic 
inflation J35| was driven by </>, we find M ~ lO^^GeV and A < 10"^'^ [36]. Then in order 
that radiative corrections do not disturb the potential we require h < 10^^, / < 10^^ and 
M < lO^^^GeV. 

Below we consider the effect of each interaction term separately. 

B. Reheating through Yukawa coupling 

First we consider the case the inflaton is coupled only with fermions ip and ip through 
Yukawa coupling. In this case preheating due to parametric resonance is unimportant and 
the dissipation rate is given from (jH^ as 



2 



f 

Tf(T) = —M 



in 



l-2np( — 



(145) 



in the perturbative regime M > f(j). If the reheat temperature turns out to be much lower 
than M the dissipation rate agrees with the conventional calculation which gives one particle 
decay rate of 0, 

rF,conv = ^M, (146) 



Stt 



which gives 



On the other hand, if M is so small that the reheating is completed in a high temperature 
regime T ^ M, we find from ()145j) that 

TphiKhiT) = . (148 

Inserting it in ()140p the reheat temperature is approximately given by 

Trm.^-[:;^J [-^^^) ■ (149) 
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This formula applies when T/j^high ^ M, or 

1/2 f2 



-«(^) 4^--°V.-(l^)-- 



(150) 



In this case relaxation of the inflaton is delayed due to Pauli blocking. But the discrepancy 
between the new result ()149|) and the conventional one ()147p is rather modest, 

with a weak dependence on the model parameters (M/~^)^/^. 

Finally we confirm consistency of our analysis. The condition /0 < M is satisfied at the 
time of reheating for 

/ f \7/2 

M > 4 X 10 V/^^ ( ^^ ) GeV. (152) 

On the other hand, the condition that thermal mass of generated by Yukawa coupling, 
/T, is smaller than M reads 

M > 5 X lOV*'^^^ ( ^^ ) ^^^- ^^^^) 

We see that ()150|) . p52|) . and ()153|) can easily be satisfied simultaneously. 

C. Reheating through three body bosonic interaction 

Next we consider the three body bosonic interaction 7V10x^, which induces a dissipation 
rate (jTTjl 



^b{T) 



M 



2 



l + 2nB[—j 



(154) 



SvrM 
In the high temperature limit the dissipation rate is enhanced as 

Tb..,.(T) = 5-gj. (155) 

This expression applies when T ^ M and M^ > M.(f). The latter requirement is the same as 
the condition broad resonance is no longer effective. We are interested in the case reheating 
is completed in this high temperature regime. Using the formula 

_ / 90 \ ^/^ I 

^ii,high — —5 — v r_B,high(^i?,high)Ms, (156) 

KT^ g*J ^ 
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we would obtain the reheat temperature 

TuM,^ - [^J ^^^^. (157) 

Consistency Tp^high ^ M would then read 

The other condition for (|155|) to apply, namely Ai(f) < M^, requires the radiation energy 
with temperature ()157j) should be smaller than M^ / M^ = P(p{(f) = M'^/Ai), which reads 

5/7 

(159) 



M>1.3xlO^V.-A^(^^) GeV. 



Clearly, (jl58|) and (jl59|) are hardly compatible with each other. 

This means that if reheating is governed by the high-temperature dissipation rate ()155|l 
the reheating process occurs shortly after the field amplitude gets smaller than M"^ / M. when 
()155|) becomes applicable. Then the use of the formula ()156|) is inappropriate and we should 
use 

90 y/^ ^^^ f?,oeV'^ M^i^ 






A<V2 
- 3,5 X lO'V.--,- (_|_y'' (To4v)"^GeV. ,160) 

where He denotes the Hubble parameter when becomes as small as M'^/Ai, namely 

Hc = ^ . (161) 

Here ,^ > 1 is a parameter which represents contribution of residual radiation energy density 
Pr created by the parametric resonance. It is defined by pr = (^^ — 1)^^ and would reduce to 
unity if preheating was totally negligible. Now the consistency condition Tr high ^ M reads 

1/2 

(162) 



M»(^) r'M = 8.igV/'r'M. 



Let us confirm the dissipation rate p55|) at the temperature p60|) is larger than He, which 
yields 

M < 1.3 X 10-,V^/^^r^/^ (l5T^) '^' ^^^' ^^63) 
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which is consistent with ()159p . because these two inequahties have been derived from the 
opposite conditions. 

On the other hand, if the conventional dissipation rate F^conv = M."^ /{9>tiM) was larger 
than Hci or 

M < 6.4 X lO-r^/^ ^^^.^^^ j GeV, (164) 

the conventional reheating process would also proceed as rapidly as to give the same reheat 
temperature. Hence the effects of the high-temperature enhancement of the dissipation rate 
is prominent only when the inequality 

'■' ' ^°"^"" (lo4v)'"°^^ < '' < '■' ^ 10'V.-'»r"' (l54v)^"GeV. 

(165) 

is satisfied. As a result the ratio of the new reheat temperature p60p to the conventional 
estimate, 

(166) 
is at most 

rfl,l,lgh_V«? M' ^„o„- 3/14,-1/7/' M ^-"" 



vJ <««•""«"" (lo4v)" ■ («^) 



T/J,conv 31/4 M^/m^J 

D. Reheating through setting-sun diagrams 

Finally we consider reheating through dissipation due to the setting-sun diagrams, in 
particular, arising from the interaction /i^0^x^/4 corresponding to the diagram L7. The 
dissipation rate which applies at high temperature T ^ M = m^ and low field amplitude 
(p < M/h after the broad resonance regime is (jl(J6|) or the second line of (I1U8|) . 

YsiT) = .21_^ (168) 

*^ ' 2567rM' ^ ^ 

which has the same temperature dependence as the Hubble parameter in the radiation 
dominated universe. Hence in order to reheat the universe completely due to this dissipation 
term we must have 
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namely, 

M < 1.9 X lOV*"'/' (1F3) '^''^- ^'^^^^ 

Let us first pretend that preheating is negligible and all the radiation conies from dissi- 
pation p68|) . Then radiation density and the temperature for (p < M/h are given by 

Here we have used flMlj) . which is not strictly valid when the dissipation rate depends on 
background temperature but still gives reasonably correct order of magnitude. Inserting 
dnH) to (HHHI) we find 

when the inequality ()170p is satisfied. Thus we find that Ts(T) is already close to H even 
if we take into account only the radiation produced by perturbative processes governed by 
p68p . and in this case it can be larger than H when 

Y^j GeV, (173) 

to reheat the universe soon after the epoch < M/h. 

If we include the effect of preheating the cosmic temperature could be higher than (jlTlj) . 
Then Ts{T) could be larger than H at the epoch (p < M/h under p70|) . For this to be 
the case, preheating is only required to create twice or more radiation than perturbative 
processes during broad resonance regime. 

Finally we examine the consistency of our analysis, /iTr < M -C Tr, where the former 
is the condition that thermal mass of generated through h"^ ((? x^ / A interaction remains 
smaller than M. Denoting the residual radiation energy density due to preheating by pr = 
(^^ — l)M'^(j? as in the previous subsection, the reheat temperature T^^s reads, 

^- - (^) '" (I) "' " - "«•"'* (i^) "'^ ^"'"' '^^^' 

because the total energy density, ptot = ^^^^0^, is efficiently converted to radiation at 
~ M/h in this scenario. Thus the desired condition is easily satisfied. 

Due to the strong dependence of the dissipation rate on the coupling constant h^, the 
above processes are operative only for inflation with a small mass scale M ~ lO^GeV for 
small coupling h ~ 10^'^. It is interesting to note, however, that in this case the scalar field 
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can dissipate its energy to get thermalized even in the absence of interactions that lead 
to decay, such as f(p4'ip or A10x^. 

VI. EVAPORATION RATE OF OSCILLATING FLAT DIRECTION 

A. Behavior of flat direction after inflation 

In generic supersymmetric theories there are a number of directions in scalar field config- 
uration space along which the potential vanishes except for a soft supersymmetry-breaking 
mass term. Such a flat direction field may acquire a large expectation value of order of Mg, 
beyond which the potential blows up exponentially in minimal supergravity, by accumu- 
lating quantum fluctuations during inflation and they start coherent field oscillation only 
after the Hubble parameter has decreased to the soft mass of order of ~TeV or so. Then 
the large-amplitude oscillation can easily violate baryon and lepton number conservation to 
generate baryon-to-entropy ratio up to 0{1). This is the original picture of the Affleck-Dine 

n n 

baryogenesis p supplemented by inflationary cosmology j37l |. 

Later the effect of finite-density supersymmetry breaking, especially, the Hubble-induced 
mass term and the importance of the nonrenormalizable terms in the superpotential were 



investigated by Dine, Randall and Thomas 38|. They included the following nonrenormal- 
izable term in the superpotential W. 

W D -^^'', (175) 

where (/? denotes a fiat direction field, n is an integer larger than 3, A„ is a constant of order 
of unity, M^ is some large cut-off scale such as the GUT or Planck scale. Together with the 
Hubble-induced mass term, the scalar potential reads 



V{ip) ~ ml\ip\'^ - cH'^\ip\ 



{Am^l2 + aH)Xnip' 



+ H.C. 



W 



2n-2 



+ |A.P^, (176) 



where A, a, and c are dimensionless quantities of order of unity, and WI3/2 ~ ITeV is the 
gravitino mass which we expect is of the same order of the soft mass m^p. It is important to 
have c negative. Then the instantaneous minimum of ip is located at 

Iv^l ^ (i7M,"-^)i/('^-2), (177) 

when H ^ m3/2 ~ Tn^. 
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In the scenario of Dine, Randall, and Thomas 38[, the scalar field starts oscillation 
with the angular frequency m^ as H becomes less than m^p ~ ""^3/25 and baryon number 
is generated. At this stage, however, there was a fear that the scalar condensate might 
evaporate before sufficient oscillation was achieved, because they postulated that the Affieck- 
Dine field would evaporate due to the scattering by thermal particles produced during the 
inflaton oscillation regime with the temperature ()142|) . As shown in [39| the scattering 
crosssection of zero-mode particle with mass m^ ~ ""^3/2 by a thermal particle such as a 
fermion with Yukawa coupling / with energy and momentum ~ T is of order of 

Pa 
^ - -^, (178) 

where a = g'^/{ATT) is a gauge coupling strength. Multiplying the number density of thermal 
particle n r^ T^, the ratio of the scattering rate F^ to H reads 



rsj 



> 



H m3/2H m3/2 
where we have neglected numerical factors, and the inequality is saturated in the radiation 
dominated regime with H ~ T'^/Mq. Apparently this quantity is much larger than unity 
for reasonable values of m3/2, /, and a. 

If the fiat direction interacts with thermal particle as above, however, its potential ac- 
quires finite-temperature corrections such as a thermal mass term ~ f'T at the same time. 
As a result the fiat direction may start coherent oscillation much earlier than previously 
assumed l40ll. Then the above estimate of the evaporation rate does not apply, and the 
authors of J40[ used the scattering rate with thermal particles, T^ ~ f^T or T^ ~ g^T, for 
the evaporation rate of the fiat direction, where g {> f) is the gauge coupling. The former 
formula applies when f\ip\ < T < g\ip\ and the latter for g\ip\ < T. 

These crude estimate has been refined by Anisimov and Dine J22|- They observed the 
center-of-mass energy between zero mode condensate with mass f'T and a massless thermal 
particle is of order f^/'^T whose square should replace the denominator of p78|) . As a result 
they find 

r^ ~ faT. (180) 

Whichever type of masses are used, in all the above estimates of the evaporation rate of 
the oscillating fiat directions, it was analyzed with a picture of particle-particle scattering. 
However, since the zero-mode field oscillation occupies the entire space homogeneously, it 
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would be more appropriate to regard it as a coherent condensate rather than particles. Hence 
we should use the formalism developed in the present paper instead. 

B. Dissipation rates of oscillating flat direction with a thermal mass 

Although flat direction fields are complex scalar fields, if the main driving force of their 
oscillation is their mass term, we can approximately regard them as a pair of independent 
real scalar fields and use our results based on finite-temperature nonequilibrium field theory 
to calculate the evaporation rate. If, on the other hand, y? had a large initial value ip ~ Mq 
and the condensate acquired huge baryon number density, its eva por ation would be delayed 
because chemical potential of bosons cannot exceed their mass (4lL l42l |. We assume that 
initial value of ip is regulated to a sufficiently small value ()177|) due to the nonrenormalizable 
terms in the potential (jl76|) and consider the situation the flat directions dissipate their 
energy through the relevant dissipation rate we have obtained in §111 and §IV. 

These fields can possess all types of interactions discussed so far, namely, Yukawa cou- 
pling fifiipip, three-body scalar interaction Aiipx^ and biquadratic interaction /i^|(^p|xp. 
Here typical value of Ad is ffi with fi being the energy scale of the standard model that 
emerges as the coefficient of HuH^ term in the superpotential of the minimal supersymmet- 
ric standard model, while we expect several types of biquadratic interactions with Yukawa 
coupling strength h = f and gauge coupling strength h = g > f. If the cosmic temperature 
is higher than g\(f\ we expect (p has a thermal mass of ~ gT and it drives coherent oscillation 
when H < gT. For f\ip\ < T < g\ip\ the flat direction has a thermal mass ~ f'T and it can 
also drive oscillation when H < fT. Here we flrst write down the dissipation rates from 
various interactions for each case and then consider which rates are applicable in the next 
subsection. 

First we consider the case T > g\(f\ and H < gT so that if is oscillating with the angular 
frequency M = gT < T. Using dHHI), (EB, dlDHl), and (IHij) . we can list the rate of each 
dissipation channel together with the range of its applicability. 

Tf = - — - = ^^^T, for fU < gT, (181) 

327rT 327r ' J iri y , k j 

for fM < g^T\ (182) 



. ^ (//i)^T _ /V 
^ 27rM2 2Tig^T' 
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^-— = ^-T, for U < T, 

U-kM 647r ' iv-l ~ ' 

~ TT^TT^^ for |y?| > T, 



(183) 



L 1287r2m^ 1287r2|y;| ~ 1287r2' 

^ TTT^' for IV^I < -^' 



M 



USiiM USiiT ~ 12877' 
fV 9f^T fV^ 



;i84) 



I. 1287rM 1287r ~ 1287r^T 

In the second equality of ()183|1 we have put m^ = g\ip\. If two or more channels are at work, 
the total dissipation rate is given by their sum. 

Next for f\{p\ < T < g\ip\ (< \{p\), fields coupled to ip with gauge coupling strength are 
not thermalized and only those coupled with Yukawa coupling strength are relevant. Hence 
when H < fT, the scalar field oscillates with the angular frequency M = fT ^ T. In this 
case only the following two channels could be nonvanishing. 

r^ ^ iM!j = JfL, for /x|<^| < fT', (185) 



fT' ^ gfT 



1287r2|v9| ~ 1287r2 
C. Dissipation rate at the onset of field oscillation 



;i86) 



Finally we combine the above results with the thermal history and the initial condition of 
(p after inflation in order to evaluate the dissipation rate at the onset of field oscillation. After 
inflation, ip is expected to trace the instantaneous minimum ()177|) \(p\ ~ (HM^~^Y^^'^~'^ 
until the onset of field oscillation due to a thermal mass. For definiteness let us consider 
the case preheating is not effective so that the cosmic temperature is given by ()143|) . T ~ 
{T'^HMgY'^, during the inflaton field oscillation regime. Let us also take n = 4 below. 

The flat direction starts oscillation with a frequency gT if both T > g\ip\ and gT > H 
hold true, or with a frequency fT when both T > f\ip\ and fT > H hold, whichever comes 
earlier. The condition T > g\ip\ is satisfied when 

H < g-'TlM;'MG = i/th,g, (187) 
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while T > f\(f\ applies when 

H < f-^TlM;^MG = //th,f. (188) 

On the other hand, the inequality gT > H holds when 

H < g^/^T^J^M'f = if,,g, (189) 

while fT > H is fulfilled when 

H < f/^Tl^^M^^^ = H,^i. (190) 

We find iJs,f < -f^s.g and i^th,g < -f^th.f- 

The fiat direction starts oscillation with the frequency M = gT at H = niin(ifs^g, Hth,g), 
if min(ifs g, ii/^th,g) > niin(ifs_f, i^th,f)- This inequality holds true if 

Tr > g^fM^J^M-'/' = 8 X 10^° f^Y ( ^) ( ^f* ) GeV. (191) 

In this case the dissipation rate is given by ()181|) through p84p depending on the value of 
\ip\ and T then. 

On the other hand, if min(ifs_g, i^th,g) < ^^{Hgf, Hth,{), the scalar field starts oscillation 
with the frequency M = fT aX H = niin(iJs,f5 -f^th.f)- This happens if 

Tr < g^fMl'^M-^'\ (192) 

and the dissipation rate is given by p85p or ()18(ij) . 

To conclude we have calculated the evaporation rate of the fiat direction at the onset 
of its oscillation for n = 4. We find that in some cases the rate may be larger than the 
previous estimate based on particle-particle scattering picture ()180|) but the time scale of 
evaporation is long enough that significant oscillation is certainly possible before evaporation. 
Once it starts oscillation, the evolution of the field amplitude \(p\ becomes different from 
()177j) . so one must solve its evolution together with the thermal history after inflation in 
order to determine when the flat direction completes thermalization. This issue will be 
analyzed elsewhere together with the amount of baryon asymmetry produced, where two- 
loop logarithmic correction to the effective potential J39|, which we have neglected here, will 
also be included. 
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VII. DISCUSSION 

In the present paper we have developed a formahsm to investigate the relaxation processes 
of an oscillating scalar field interacting with various particles in a thermal state using 
the in-in formalism of nonequilibrium quantum field theory. Integrating out those thermal 
particles interacting with 0, we have obtained the effective action for which is complex 
even if it is a real scalar field. This is a result of coarse-graining and manifestation of the 
dissipative effect on to those integrated out. The real equation of motion is obtained by 
introducing auxiliary fields, C,a{x) and ^m{x), which act as an additive and a multiplicative 
noise term, respectively. The former originates from interactions linear in such as Yukawa 
coupling or three-body bosonic interaction, while the latter is from quadratic or higher-order 
interactions in 0. It induces noises on the effective mass of 0. 

The equation of motion has terms nonlocal in both space and time as a result of quantum 
corrections. In the linear regime when higher-order terms in are negligible in the equation 
of motion, these nonlocalities can easily be handled because its Fourier modes are decoupled 
from each other. As a result we can find an analytic solution for each mode from which we 
can extract the dissipation rate. On the other hand, the dissipation rates from multiplicative 
interactions are read from the equation of motion itself. 

Quite generally, the memory kernels, which generate nonlocal terms in the equation 
of motion, are determined by the imaginary part of the Green functions relevant to each 
diagram, while the noise correlation functions are identical to the real part of the same 
function up to a numerical factor. We have found that for all the interactions discussed 
here, the Fourier transform of the memory kernel and that of noise correlation function take 
a specific ratio which is determined only by the temperature and the angular frequency of the 
mode. This relation is achieved by microphysical detailed balance relation. It also leads to 
the well-known fiuctuation-dissipation theorem for low-momentum modes, which guarantees 
that the scalar field relaxes to a state the equipartition law is satisfied. For higher- momentum 
modes the scalar field relaxes to the thermal equilibrium state with the same temperature 
where the number density of each quanta consists of the boson distribution function and its 
zero-point vacuum component. 

Although we have shown the fate of the oscillating scalar field is the same equilibrium 
state, the time scale of relaxation to it is strikingly different depending on the nature of 
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interactions. In the case of Yukawa coupling with fermions, the dissipation rate takes a 
smaller value at finite temperature than the zero-temperature decay width due to the Pauli 
blocking. On the other hand, in the case of bosonic three-body interaction, the dissipation 
rate is larger than the zero temperature decay rate due to the induced effect. As a result we 
have seen the reheat temperature after inflation may be somewhat changed from conventional 
estimates, and that in an extreme case the inflaton can dissipate its energy even without 
linear interactions that leads to its decay. 

The temperature dependence on the dissipation rate may also affect on the property and 
the spectrum of density fluctuations. It has been known for a long time that primordially 
isocurvature fluctuations that were stored in a long-lived scalar fleld during inflation when 
it was subdominant are converted to the adiabatic ones as its energy density tends to dom- 
inate the Universe late r | 43| . Such a property has been utilized in some models of non-scale 
invariant fluctuations (44l . l4q . Nowadays the above conversion mechanism from isocurva- 
ture to adiabatic fluctuations is called the curvaton scenario |46|. When the curvaton fleld 
decays, the dependence of their dissipation rate on the background temperature may in- 
duce additional fluctuations just as in the modulated coupling scenario [47]. On the other 
hand, a model of baryogenesis has been proposed in which small fluctuation in the inflaton's 
dissipation rate induces enhanced baryon-number fluctuations '48]. These possible effects 
on density fluctuations due to the temperature dependence on the dissipation rate will be 
studied elsewhere. 

The dissipation associated with interactions linear in such as Yukawa coupling and 
three-body bosonic coupling can be interpreted in terms of decay, while the setting-sun 
diagram from biquadratic coupling K^cj^x^ /A and quartic coupling A0'^/4! induces dissipation 
associated with scattering. The former is suppressed when the would-be decay products 
are more massive than the oscillation frequency, but the latter is effective even in this 
regime. Consequently the dissipation rate of the Affleck-Dine flat direction fleld shows a 
rather complicated behavior depending on the evolution of its oscillation amplitude and the 
temperature. Although the dissipation time scale is much longer than the oscillation period, 
whether sufficient baryon number is generated or not depends on the magnitude^ of the A 
terms as well, in particular on the presence or the absence of the thermal A term |39|. Hence 
the both ingredients should be analyzed properly to yield final baryon asymmetry. 

In the present paper we have concentrated on the fate of the zero-mode oscillation but 
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we can easily obtain the dissipation rates of higher-momentum modes using our formula and 
we expect they have larger dissipation rate. This may affect formation of Q-balls. 

Thus there are a number of interesting problems associated with dissipation of fiat direc- 
tions remaining. We hope to return to these issues in near future. 
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